Measuring neutron-star properties via gravitational waves from neutron-star mergers 
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We demonstrate by a large set of merger simulations for symmetric binary neutron stars (NSs) 
that there is a tight correlation between the frequency peak of the postmerger gravitational-wave 
(GW) emission and the physical properties of the nuclear equation of state (EoS), e.g. expressed 
by the radius of the maximum-mass Tolman-Oppenheimer-Volkhoff configuration. Therefore, a 
single measurement of the peak frequency of the postmerger GW signal will constrain the NS EoS 
significantly. For optimistic merger-rate estimates a corresponding detection with Advanced LIGO 
is expected to happen within an operation time of roughly a year. 
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The properties of high-density matter as in the cores of 
NSs, in particular the EoS, are still incompletely known, 
because the physical conditions are not directly accessi- 
ble by experiments. Theoretical models for supernuclear 
matter are ambiguous and suffer from uncertainties of 
nuclear data required as input for these calculations [l[. 

NS properties are intimately linked to the adopted EoS 
because the latter determines the stellar structure by the 
Tolman-Oppenheimcr- Volkhoff (TOV) equations 0, Q. 
Hence, constraints on the NS EoS can be deduced from 
astrophysical observations (e.g. 0]), as alternatives to nu- 
clear models 0] and laboratory experiments [H . 

NS mergers may also yield information about the nu- 
clear EoS, because the dynamics of the coalescence de- 
pend sensitively on the behavior of high-density matter 
(see 0, Q for reviews). Consequently, the EoS leaves 
an imprint on the GW signal of NS mergers. However, 
the systematic dependences of the inverse problem, i.e. 
which EoS (or NS) properties can be derived from a par- 
ticular GW detection, are still not completely explored 
(see 0,0413] and refs. therein). In this letter we report 
on a tight correlation between NS parameters and thus 
EoS characteristics and the dominant frequency of the 
postmerger GW emission revealed by a systematic study 
with 18 microphysical EoSs. Our survey is in particular 
important because the second- gene ration interferometric 
GW detectors of Adv. LIGO [l7| and Adv. Virgo [3 
go into operation within the next years. NS binaries are 
considered a major target of these instruments with an 
estimated detection rate of 0.4 to 400/yr (l9j . 

Our simulations are performed with a 3-D relativis- 
tic smoothed particle hydrodynamics (SPH) code, which 
solves the Einstein field equations assuming conformal 
flatness and employing a GW backreaction scheme within 
a post- Newtonian framework [45[ [2(| 21]. The imple- 



mentation allows the usage of tabulated microphysical 
EoSs including thermal effects, or arbitrary barotropic 
EoSs (e.g. zero-temperature EoSs for equilibrium to weak 
interactions, the so-called /3-cquilibrium). The latter are 
supplemented by an ideal-gas component with an ideal- 
gas index T t h = 2 to mimic thermal effects j22[. 



The calculations start from quasi-equilibrium orbits 
about two revolutions before the merging of the NSs, 
which are assumed to be initially cold and in neutrino- 
less /3-equilibrium. Because tidally locked binaries are 
unlikely to occur [23j, the stars are set up as nonrotat- 
ing, which is a valid approximation even for millisecond 
NSs, whose rotation is still slow compared to the orbital 
period. If not noted otherwise the NSs are modeled by 
about 340,000 SPH particles. 

In total we employ 18 different microphysical EoSs (see 
Tab. U for the nomenclature and references). Seven of 
these EoSs include thermal effects consistently. The re- 
maining ones describe nuclear matter at zero tempera- 
ture and are labeled with "+Pth" in Tab. |U The mass- 
radius (M-R) relations, the maximum masses M max of 
nonrotating NSs and the corresponding (minimum) radii, 
denoted as i? max , for all used EoSs are shown in Fig. Q] 
The maximum-mass configurations (Tab. UJ are marked 
by symbols. The scatter in Fig. [T] illustrates the diversity 
of the microphysical models underlying our study. 

We consider EoSs with M max in the range of 1.80 Mq 
to 2.76 M and i? max from 9.30 km to 14.30 km with- 
out any special selection procedure except that we re- 
quire A/ max > 1.8 M Q . The lower limit of 1.8 M is 
motivated by the detection of a pulsar with a mass of 
(1.97 ± 0.04) M Q 0. Although this observation rules 
out some EoSs of our sample, we do not disregard these 
models, because at lower densities (as present in 1.35 Mq 
NSs and in the merger remnant where strong rotational 
and thermal effects come into play) these EoSs may still 
provide a viable description of nuclear matter. Further- 
more, the inclusion of these EoSs demonstrates the valid- 
ity of the relations between merger and EoS properties 
discussed below over a wider parameter range. 

For each EoS listed in Tab.[D we simulate the merger of 
two stars with 1.35 Mq. This setup is chosen because pul- 
sar observations and population synthesis studies suggest 
these systems to be most abundant [39|. After energy 
and angular momentum losses by GWs have driven the 
inspiral of the NSs for several 100 Myrs, there are two dif- 
ferent outcomes of the coalescence. Either the two stars 
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TABLE I: Used EoSs. M max and i? max are mass and radius 
of the maximum-mass TOV configuration, / pca k is the peak 
frequency of the postmerger GW emission with the FWHM (a 
cross indicates prompt collapse of the remnant). //i z (/ pea k) 
is the effective peak amplitude of the GW signal at a polar 
distance of 20 Mpc. The tables of the first five and next seven 
EoSs are taken from [24| and [25|, respectively. 



EoS with 
references 



Sly4 [26] +r th 
APR [27] +r th 
FPS [28] +r th 
BBB2 [29] +r th 
Glendnh3 [30j+r th 
eosAU [31] +r th 
eosC [32] +r th 
eosL [33] +r th 
eosO [34] +r th 
eosUU [31] +r th 
eosWS [31] +r th 
SKA [35] +r th 
Shen [36] 
LS180 [35] 
LS220 [35] 
LS375 [35] 
GS1 [37] 
GS2 [381 



M max 
[M ] 



2.05 
2.19 
1.80 
1.92 
1.96 
2.14 
1.87 
2.76 
2.39 
2.21 
1.85 
2.21 
2.24 
1.83 
2.04 
2.71 
2.75 
2.09 



R max 

[km] 



10.01 

9.90 

9.30 

9.55 

11.48 

9.45 

9.89 

14.30 

11.56 

9.84 

9.58 

11.17 

12.63 

10.04 

10.61 

12.34 

13.27 

11.78 



/ poak , FWHM 
[kHz] 



3.32 
3.46 

x 

3.73 
2.33 

x 

3.33 
1.84 
2.66 
3.50 

x 

2.64 
2.19 
3.26 
2.89 
2.40 
2.10 
2.53 



0.20 
0.18 

0.22 
0.13 

0.22 
0.10 
0.11 
0.17 

0.13 
0.15 
0.25 
0.21 
0.13 
0.12 
0.12 



/M/p 

oak , 

[1Q- 21 ] 



2.33 
2.45 
x 

1.33 
1.27 

x 

1.27 
1.38 
2.30 
2.64 
x 

1.96 
1.43 
1.19 
1.63 
1.82 
1.46 
2.15 



directly form a black hole (BH) shortly after they fuse 
("prompt collapse"), or the merging leads to the forma- 
tion of a differentially rotating object (DRO) that is sta- 
bilized against the gravitational collapse by rotation and 
thermal pressure contributions. Continuous loss of angu- 
lar momentum by GWs and redistribution to the outer 
merger remnant will finally lead to a "delayed collapse" 
on timescales of typically several 10-100 ms depending 
on the mass and the EoS. For EoSs with a sufficiently 
high M max stable or very long-lived rigidly rotating NSs 
are the final product. 

A prompt collapse occurs for three EoSs of our sample 
(marked by x in Tab. Q] and Fig. [T]). One observes this 
scenario only for EoSs with small R max . In the simula- 
tions with the remaining EoSs DROs are formed. The 
evolution of these mergers is qualitatively similar. The 
dynamics are described in [HED]. 

For all models that produce a DRO the GW signal is 
analyzed by a post-Newtonian quadrupole formula poj ]. 
The inset of Fig. [2] shows the GW amplitude of the plus 
polarization at a polar distance of 20 Mpc for NSs de- 
scribed by the Shen EoS. Clearly visible is the inspi- 
ral phase with an increasing amplitude and frequency 
(until 5 ms), followed by the merging and the ringdown 
of the postmerger remnant (from 6 ms). All DROs are 




12 14 
R [km] 

FIG. 1: NS M-R relations for all considered EoSs. Red curves 
correspond to EoSs that include thermal effects consistently, 
black lines indicate EoSs supplemented with a thermal ideal 
gas. The horizontal line corresponds to the 1.97 Mq NS @]. 



stable against collapse well beyond the complete damp- 
ing of the postmerger oscillations. In Fig. [5] we plot 
the spectra of the angle-averaged effective amplitude, 
/i av = 0Afh z (f) (see e.g. 0), at a distance of 20 Mpc 
for the Shen EoS (solid black) and the eosUU (dash- 
dottedVtogether with the anticipated sensitivity for Adv. 
LIGO [H and the planned Einstein Telescope (ET) 0. 



Here h z (f) = y {\h + \ 2 + \h x | 2 )/2 is given by the Fourier 

transforms, h + / y , of the waveforms for both polariza- 
tions observed along the pole. As a characteristic feature 
of the spectra a pronounced peak at / poa k = 2.19 kHz for 
the Shen EoS and 3.50 kHz for eosUU is found, which is 
known to be connected to the GW emission of the merger 
remnant Q. Recently, this peak has been identified as 
the frequency of the fundamental quadrupolar fluid mode 
(f-mode) [41| . For all models producing a DRO the spec- 
tra are sharply peaked in the kHz range around / pC ak with 
a FWHM below 250 Hz. Values of / pca k, the FWHM and 
the height of the peak for all models are listed in Tab. fl] 
For the Shen EoS Fig. [2] also shows the results of a run 
starting 3.5 revolutions before merging (red line), for a 
calculation with 1,270,000 SPH particles (blue), and for 
a simulation neglecting the GW backreaction in the post- 
merger phase (green) confirming the insensitivity to these 
choices. The initial rotation state of the NSs is known 
to affect /peak only insignificantly 12j. Furthermore, our 



/peak values agree within a few per cent with the results 
of fully relativistic simulations (e.g. 3.35 kHz for the 
APR EoS in [Til] ) - The uncertainties associated with the 
r t h— ansatz for thermal effects are below 10 per cent [22| . 

Our systematic study reveals that the peak frequency 
/peak of the GW signal produced by the oscillating, hot, 
highly deformed DRO is determined by characteristic 
properties of NSs on the AI-R- sequence for nonrotat- 
ing TOV solutions. In Fig. [3] / pca k is plotted against 
-Rmax (crosses and triangles) and an obvious empirical 
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FIG. 2: Orientation-averaged spectra of the GW signal for 
the Shen (solid) and the eosUU (black dashed-dotted) EoSs 
and the Adv LIGO (red dashed) and ET (black dashed) unity 
SNR sensitivities. The inset shows the GW amplitude with + 
polarization at a polar distance of 20 Mpc for the Shen EoS. 



correlation is visible. / pea k is higher for smaller i? max - 
The outlier (triangle) belongs to the simulation for the 
Glendnh3 EoS, which has a strikingly different M-R re- 
lation (dashed line in Fig. [1]), which seems in conflict 
with theoretical knowledge of EoS properties at subnu- 
clear densities Ignoring the outlier, the remaining 
"accepted models" exhibit an even stronger / pe ak--Rma X 
correlation (line in Fig. [3]) . Already one determination 
of /peak could therefore seriously constrain the M-R re- 
lation and consequently the nuclear EoS. Additionally, 
simulated mergers of 1.2 M© -1.5 M© binaries for se- 
lected EoSs (circles) demonstrate that the relation be- 
tween /p ea k and i? m ax is not very sensitive to the ini- 
tial mass ratio [12]. Squares in Fig. [3] display results 
for 1.2 Mg-1.2 Mq mergers. For those / p0 ak is clearly 



lower [12| with differences being larger for smaller R max . 
But also for the symmetric binaries with lower mass a cor- 
relation seems to exist. We stress that the total binary 
mass M to t is measurable by the GW inspiral signal (42j . 

/peak turns out to correlate also with other proper- 
ties of static, cold NSs: From Fig. 2] (left panel) a close 
relation between the radius R1.35 of a 1.35 M© star 
and /peak is evident. A similar coupling is found be- 
tween /peak and the maximum central density p max of 
nonrotating NSs, where higher p max yield higher /peak- 
However, no clear correlation exists between / pca k and 
M max , though typically a lower M max gives a higher 
/peak, and /peak > 2.8 kHz seems incompatible with 
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A tight relation between the frequency of nonlinear 
oscillations of the hot, highly deformed DRO and the 
properties of cold, static, spherical NSs (on the M-R- 
sequence) is an empirical finding of this work. This links 
to numerical calculations which showed that for any EoS 
the frequency of the f-mode (generating the GW radia- 
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FIG. 3: Peak frequency of the postmerger GW emission vs. 
radius of the maximum-mass TOV solution. Blue cases are 
excluded by 0] . Error estimates are based on a Fisher matrix 
analysis for a source at 20 Mpc. The line is a least square 
fit -Rmax oc f~ci^ f° r t ne accepted models. The triangle is an 
outlier (see text) and the squares correspond to models with 
lower Mtot- See text for other symbols. 
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FIG. 4: Peak frequency vs. radius of a 1.35 M e NS (left) and 
vs. -J A/tot / -Rmax in geometrical units (right) with M to t being 
the binary mass. Symbols have same meaning as in Fig. 



tion at /peak 14l|) depends nearly linearly on the square 
root of the mean density, (M/R 3 ) 1 / 2 , reflecting the dy- 
namical time scale of a stellar object [43J . While M here 
is identified with the mass of the DRO, approximately 
given by M to t, the empirical correlation of Fig. [3l how- 
ever, means that R in this formula can be expressed by 
Rmax of the maximum-mass TOV configuration. With 
Mot being fixed, this means that / pea k oc -R max 5 , which 
is verified by the right panel of Fig. @J where except for 
the mentioned outlier a clear linear scaling is visible. A 
fit to -Rmax (/peak) oc (line in Fig. [3j, using only 

accepted models, reveals a maximum residual of 300 m. 
The tight correlation of / pca k and (M to t/i? max ) 1/2 im- 
plies that the radius of the DRO and i? max are strongly 
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linked. Such a strong link has indeed been empirically 
found to exist between i? max and the radius of the most 
massive, rigidly rotating NS for any EoS [44 |, and seems 
to exist also for differentially rotating NSs with 2.7 M . 

Using the postmerger signal alone and correcting the 
underestimation of 40 % of the GW amplitude by the 
quadrupole formula [lo| , a SNR of 2 (the inclusion of the 
inspiral signal increases the SNR significantly) yields an 
optimal detection horizon of about 20-45 Mpc (depen- 
dent on the EoS) for Adv. LIGO. This corresponds to 
145-1190 Milky Way Equivalent Galaxies accessible for a 
GW search, taking into account the reduction due to ran- 
dom source location and orientation [l!| . The "realistic" 
and the "high" merger rates of [l9| imply a detection rate 
of 0.01-1 events/yr for Adv. LIGO. With the planned 
ET [4(J and its higher sensitivity several observations of 
/peak per year will become very likely. 

For polar distances of 20 Mpc / poa k-uncertainties of 
typically 50 Hz and at most 160 Hz can be estimated 
from the Fisher information matrix for neighboring cases 
of accepted models following [l3j |. Corresponding uncer- 
tainties are indicated in Fig. [3] by averages for contiguous 
pairs of models. Considering in addition the residuals to 
the fits of the relations of Figs. [3]and|4l a measurement 
of /peak will allow to determine the NS radius with an 
accuracy of several 100 m. These prospects are compa- 
rable with the 1 km accuracy of the radius estimation for 
the initial NSs from the inspiral GW signal of symmet- 
ric binaries suggested in [l3| for events within a maximal 
distance of 20-100 Mpc. Both will set strong constraints 
on the M-R relation and thus the EoS @] ■ Our approach, 
however, is an independent, complementary one, also to 
the possibility of gaining EoS information from the weak 
correlation between M max and the threshold total binary 
mass Mthres that distinguishes prompt (M to t > -Mthres) 
from delayed (M to t < M t hres) BH formation [l6|, whose 
determination requires more than one GW detection Q . 

Future numerical studies should vary M to t and con- 
firm our findings by more sophisticated models of bi- 
nary mergers, e.g. considering magnetic fields, neutrino 
physics, and full general relativity Also the detectabil- 
ity of /peak should be explored in more detail, e.g. by a 
detector network. Finally, our explanation should be ex- 
amined more closely to develop a precise understanding 
of the presented correlations. 
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